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Abstract 

The "in-in" formalism is reviewed and extended, and applied to the calcula- 
tion of higher-order Gaussian and non-Gaussian correlations in cosmology. 
Previous calculations of these correlations amounted to the evaluation of 
tree graphs in the in-in formalism; here we also consider loop graphs. It 
turns out that for some though not all theories, the contributions of loop 
graphs as well as tree graphs depend only on the behavior of the inflaton 
potential near the time of horizon exit. A sample one-loop calculation is 
presented. 
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I. INTRODUCTION 



The departures from cosmological homogeneity and isotropy observed 
in the cosmic microwave background and large scale structure are small, 
so it is natural that they should be dominated by a Gaussian probability 
distribution, with bilinear averages given by the terms in the Lagrangian 
that are quadratic in perturbations. Nevertheless, there is growing interest 
in the possibility of observing non-Gaussian terms in various correlation 
functions,^ such as an expectation value of a product of three temperature 
fluctuations. It is also important to understand the higher-order corrections 
to bilinear correlation functions, which appear in Gaussian correlations. 

Until now, higher-order cosmological correlations have been calculated 
by solving the classical field equations beyond the linear approximation. 
As will be shown in the Appendix, this is equivalent to calculating sums 
of tree graphs, though in a formalism different from the familiar Feynman 
graph formalism. For instance, Maldacena^ has calculated the non-Gaussian 
average of a product of three scalar and/or gravitational fields to first order 
in their interactions, which amounts to calculating a tree graph consisting 
of a single vertex with 3 attached gravitational and/or scalar field lines. 

This paper will discuss how calculations of cosmological correlations can 
be carried to arbitrary orders of perturbation theory, including the quantum 
effects represented by loop graphs. So far, loop corrections to correlation 
functions appear to be much too small ever to be observed. The present 
work is motivated by the opinion that we ought to understand what our 
theories entail, even where in practice its predictions cannot be verified 
experimentally, just as field theorists in the 1940s and 1950s took pains to 
understand quantum electrodynamics to all orders of perturbation theory, 
even though it was only possible to verify results in the first few orders. 

There is a particular question that will concern us. In the familiar calcu- 
lations of lowest-order Gaussian correlations, and also in Maldacena's tree- 
graph calculation of non-Gaussian correlations, the results depended only 
on the behavior of the unperturbed inflaton field near the time of horizon 
exit. Is the same true for loop graphs? If so, it will be possible to calcu- 
lated the loop contributions with some confidence, but we can learn little 
new from such calculations. On the other hand, if the contribution of loop 
graphs depends on the whole history of the unperturbed infiaton field, then 
calculations become much more difficult, but potentially more revealing. In 
this case, it might even be that the loop contributions are much larger than 
otherwise expected. 
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The appropriate formalism for dealing with this sort of problem is the 
"in-in" formalism originally due to Schwinger.^ Schwinger's presentation is 
somewhat opaque, so this formalism is outlined (and extended) in an Ap- 
pendix. In section II we summarize those aspects of this formalism that 
are needed for our present purposes. Section III introduces a class of the- 
ories to serve as a basis of discussion, with a single inflaton field, plus any 
number of additional massless scalar fields with only gravitational interac- 
tions and vanishing expectation values. In Section IV we prove a general 
theorem about the late time behavior of cosmological correlations at fixed 
internal as well as external wave numbers. Section V introduces a class 
of unrealistic theories to illustrate the problems raised by the integration 
over internal wave numbers, and how these problems may be circumvented. 
In Section VI we return to the theories introduced in Section III, and we 
show that the conditions of the theorem proved in Section IV are satisfied 
for these theories. This means that, to all orders of perturbation theory, if 
ultraviolet divergences cancel in the integrals over internal wave numbers, 
then cosmological correlations do indeed depend only on the behavior of the 
unperturbed inflation field near the time of horizon exit in the cases studied. 
We can also find other theories in which this result does not apply, as for 
instance by giving the additional scalar fields a self-interaction. Section VII 
presents a sample one-loop calculation of a cosmological correlation. 

II. THE "IN-IN" FORMALISM IN COSMOLOGY 

The problem of calculating cosmological correlation functions differs 
from the more familiar problems encountered in quantum field theory in 
at least three respects: 

• We are not interested here in the calculation of S-matrix elements, but 
rather in evaluating expectation values of products of fields at a fixed 
time. 

• Conditions are not imposed on the fields at both very early and very 
late times, as in the calculation of S-matrix elements, but only at 
very early times, when the wavelength is deep inside the horizon and 
according to the Equivalence Principle the interaction picture fields 
should have the same form (when expressed in terms of metric rather 
than co-moving coordinates) as in Minkowski spacetime. 

• Although the Hamiltonian H that generates the time dependence of 
the various quantum fields is constant in time, the time-dependence 
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of the fluctuations in these fields are governed by a fluctuation Hamil- 
tonian H with an explicit time dependence, which as shown in the 
Appendix is constructed by expanding H around the unperturbed so- 
lution of the field equation, and discarding the terms of first order in 
the perturbations to the fields and their canonical conjugates. 

Given a fluctuation Hamiltonian H, we want to use it to calculate expec- 
tation values of some product Q{t) of field operators, all at the same time t 
but generally with different space arguments. As discussed in the Appendix, 
the prescription of the "in-in" formalism is that 



{Qit)) 



Texp i 



Hi{t)dt 



Q\t) 



Texp 



Hi{t)dt 



(1) 

Here T denotes a time-ordered product; T is an anti-time-ordered product; 

is the product Q in the interaction picture (with time-dependence gen- 
erated by the part of H that is quadratic in fluctuations); and Hj is the 
interaction part of H in the interaction picture. (This result is different 
from that originally given by Maldacena^ and other authors^, who left out 
the time-ordering and anti-time-ordering, perhaps through a typographical 
error. However, this makes no difference to first order in the interaction, 
which is the approximation used by these authors in their calculations.) We 
are here taking the time to at which the fluctuations are supposed to behave 
like free fields as to = — oo, which is appropriate for cosmology because at 
very early times the fluctuation wavelengths are deep inside the horizon. 

Eq. (1) leads to a fairly complicated diagrammatic formalism, described 
in the Appendix. Unfortunately this formalism obscures crucial cancella- 
tions that occur between different diagrams. For our present purposes, it is 
more convenient to use a formula equivalent to Eq. (1): 



{Qit)) = E / dtN dtN-1 ■■■ dti 

nr J —OO J — OO J —OO 

[Hi{ti)\Hi{t2),- ■ ■[Hi{tN),Q\t) 



N=0 
X 



(2) 



(with the N = term understood to be just {Q^{t))). This can easily be 
derived from Eq. (1) by mathematical induction. Obviously Eqs. (1) and 
(2) give the same results to zeroth and first order in Hj. If we assume that 
the right-hand sides of Eqs. (1) and (2) are equal for arbitrary operators Q 
up to order N — 1 in Hj, then by differentiating these equations we easily 
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see that the time derivatives of the right-hand sides are equal up to order 
N. Eqs. (1) and (2) also give the same results for t — — oo to all orders, so 
they give the same results for arbitrary t to order N. 

III. THEORIES OF INFLATION 

To make our discussion concrete, in this section we will take up a partic- 
ular class of theories of inflation. The reader who prefers to avoid details of 
specific theories can skip this section, and go on immediately to the general 
analysis of late-time behavior in the following section. 

In this section we will consider theories of inflation with two kinds of 
matter fields : a real scalar field (/^(x, t) with a non-zero homogeneous ex- 
pectation value (p{t) that rolls down a potential V{ip), and any number 
of real massless scalar fields cJn(x, t), which have only minimal gravitational 
interactions, and are prevented by unbroken symmetries from acquiring vac- 
uum expectation values. The real field ip serves as an infiaton whose energy 
density drives inflation, while the cr„ are a stand-in for the large number of 
species of matter fields that will dominate the effects of loop graphs on the 
correlations of the infiaton field.** 

We follow Maldacena,^ adopting a gauge in which there are no fluctua- 
tions in the infiaton field, so that </'(x, t) = (p{t), and in which the spatial 
part of the metric takes the form*** 

gij = a^e^'^[exp-/]ij, 7^ = 0, daij = . (3) 

where a(t) is the Robertson- Walker scale factor, 7jj(x, i) is a gravitational 
wave amplitude, and Ci'^i't) is a scalar whose characteristic feature is that 

** Standard counting arguments show that in these theories the number of factors of 8wG 

in any graph equals the number of loops of any kind, plus a fixed number that depends 
only on which correlation function is being calculated. Matter loops are numerically 
more important than loops containing graviton or infiaton lines, because they carry an 
additional factor equal to the number of types of matter fields. 

***! am adopting Maldacena's notation, but the quantity he calls ^ is more usually called 
TZ. To first order in fields, the quantity usually called is defined as — — H5p/p, while 
the quantity usually called TZ is defined as — vl* + HSu. (Here the contribution of scalar 
modes to gtj is written in general gauges as — 2a^(^Jy + d^'^' /dx^dx-'), while Sp and p 
are the perturbation to the total energy density and its unperturbed value, while 5u is 
the perturbed velocity potential, which for a single infiaton field \s 5u = —5if/(p.) In the 
gauge used by Maldacena and in the present paper Su = ^' = 0, so since C, is defined here 
as vE" to first order in fields, it corresponds to the quantity usually called TL. Outside the 
horizon TZ and C, are the same. 
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it is conserved outside the horizon,^ that is, for physical wave numbers that 
are small compared with the expansion rate. The same is true of 7ij . 

The other components of the metric are given in the Arnowitt-Deser— 
Misner (ADM) formalism^ by 
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-N^ + gijN'N^, 9io = 9ijN^, 



(4) 



where N and iV* are auxiliary fields, whose time-derivatives do not appear 
in the action. The Lagrangian density in this gauge (with SttG = 1) is 



^ 2 



-1 ^2 



where 



+ E - ^'^^<ynf - Aa-2e-2^[exp (-7)]*^' ^ diGndjU^ 

n 

Eij = l{9ij-ViNj-VjNi) , 



(5) 



(6) 



and bars denote unperturbed quantities. All spatial indices i, j, etc. are 

lowered and raised with the matrix gij and its reciprocal; Vj is the three- 
dimensional covariant derivative calculated with this three-metric; and R^^^ 
is the curvature scalar calculated with this three-metric: 



-2 -2C 

a e ^ 



The auxiliary fields N and N'^ are to be found by requiring that the action 
is stationary in these variables. This gives the constraint equations: 



n 



(7) 



Ri3) _2V- a-2e-2f [exp(-7)]*J' ^ diandjan 

n 



E^.Eh-(E\)' 



(8) 



For instance, to first order in fields (including field derivatives) the auxiliary 
fields are the same as in the case of no additional matter fields^ 



N=l + C,/H, N' 



diC + ediV-^C, 



(9) 
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where 



The fields in the interaction picture satisfy free-field equations. For ( we 
have the Mukhanov equation:^ 



+ 



Jt 



The field equation for gravitational waves is 



+ 3H 



at 



and for the matter fields 



+ - a-^V^an = . 



The fields in the interaction picture are then 

C(x,i) = / d'q [e*i-a(q)C5(i) + e-*i-a*(q)C(t) 



(11) 



(12) 



(13) 



(14) 



7i,(x,t) = / d^Yl [e*'^-"e,,(g,A)a(q,A)75(t) + e-*'i-^e*,.(g,A)a*(q,A)7*(t) 

(15) 

a„(x, t) = J d^q [e*'i-"a(q, n)a,(i) + e-'^-^a*{ci, n)a*(t)] , (16) 

where A = ±2 is a hclicity index and eij{q, A) is a polarization tensor, while 
a(q), Q!(q, A), and Q!(q, n) are conventionally normalized annihilation oper- 
ators, satisfying the usual commutation relations 



«(q) , a* (q')] = (q - q') , [«(q) , a(q') 



= . 



(17) 



a(q,A),a*(q',A') 
and 



Sxx'S''(q - q') , [a(q. A) , a(q', A')] = , (18) 



a(q,n) , Q;*(q',n') = (5nn"^^ (q - q') , a(q, ") , "(q', ^t-') =0, (19) 

Also, C<j(^)) 7ij(0) s-iid <^g(*) suitably normalized positive- frequency so- 
lutions of Eqs. (11)-(13), with replaced with — g^. They satisfy initial 
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conditions, designed to make —C,ip/H, 7jj/vT67rG, and Gn behave like con- 
ventionally normalized free fields at t — ^ — oo:^ 



H{t) ^16^ 



(27r)3/2^a(t) 



1 




(20) 



IV. LATE TIME BEHAVIOR 



The question to be addressed in this section is whether the time integrals 
in Eqs. (1) and (2) are dominated by times near horizon exit for general 
graphs. This question is more complicated for loop graphs than for tree 
graphs, such as that considered by Maldacena, because for loops there are 
two different kinds of wave number: the fixed wave numbers q associated 
with external lines, and the internal wave numbers p circulating in loops, 
over which we must integrate. It is only if the integrals over internal wave 
numbers p are dominated by values of order p ^ q that we can speak of 
a definite time of horizon exit, when q/a p/a H. In this section we 
will integrate first over the time arguments in Eq. (2), holding the internal 
wave numbers at fixed values, and return at the end of this section to the 
problems raised by the necessity of then integrating over the ps. 

There is never any problem with the convergence of the time integrals 
at very early times; all fluctuations oscillate very rapidly for q/a ^ H and 
p/a » H, suppressing the contribution of early times to the time integrals in 
Eq. (2). To see what happens for late times, when q/a <^ H and p/a ^ H, 
we need to count the powers of a in the contribution of late times in general 
loop as well as tree graphs. 

For this purpose, we need to consider the behavior of the coefficient 
functions appearing in the Fourier decompositions (14)-(16) of the fields in 
the interaction picture. In order to implement dimensional regularization, 
we will consider these coefficient functions in 2v space dimensions, returning 
later to the limit 2v 3. The coefficient functions then obey differential 
equations obtained by replacing the space dimensionality 3 in Eqs. (11)-(13) 

'^In Newtonian gauge the quantity — C(x, t)(p{t) / H (t) approaches the inflaton field fluc- 
tuation 5ip{t) for t — oo. 
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with 2z/, as well as replacing the Laplacian with —q^: 



d%{t) 



+ 



din (/),-(/) 
di 



('Kail) <f 



+ -^^m—^ + ^7,w - , 



(21) 



(22) 



(23) 



dt"^ ' dt 
At late times, when q/a ^ H, the solutions can be written as asymptotic 
expansions in inverse powers of a:^^ 



+ Cg 



t a-{t')e{t') 



dt' rt' 



roo r^i' ft too 

d.t' 



dt" 



^''{t"')e{t"') 



(24) 



dt' 2 



a2'^(t') 



a2'^(t') 



a(t") dt" 



dt" 



t" a^'^it'"' 



+ 



(25) 



a,{t) ^ a° 

/■o 

+ Eg 



,2 



dt' 



a2''(t 



— CO 

°° dt' 



a2^(t') 



a^^~\t")dt" 

^ ' Jt" a2^(t" 



(26) 



where 7^, and cr° arc the limiting values of Cqi^), lq{t)i and (Tq{t) (the 
"o" superscript stands for "outside the horizon") and C^, Dq, and Eq are 
additional constants. In any kind of inflation with sufficient expansion, the 
Robertson- Walker scale factor a will grow much faster than H or e can 



^^By t = 00 in the limits of those integrals and elsewhere in this paper, wc mean a time 
still during inflation, but sufficiently late so that a{t) is many e-foldings larger than its 
value when q/a falls below H. 
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change, and Eqs. (24)-(26) thus show that (at least for 2u > 2) the time 
derivatives of Cq, jq, and aq all vanish for q/a <^ H like 1/a^. 

If an interaction involves enough factors of ^, 7ij, and/or &„, so that 
these 1/a^ factors and any 1 / o,^ factors from the contraction of space indices 
more than compensate for the a^'^ factor in the interaction from the square 
root of the metric determinant, then the integral over the associated time 
coordinate will converge exponentially fast at late times as well as at early 
times, and therefore may be expected to be dominated by the era in which 
the wavelength leaves the horizon. For instance, the extension of Eq. (5) to 
2u space dimensions gives the interaction between a ( field and a pair of a 
fields 

jl^Caa = 7^C Y] dian diUn TTET^ X] ^i^ri diCTn 

2 n n 

^ ' n 

-fT7cE<^^ + ^cE-^ (27) 

(The Cfjcr interaction Hamiltonian given by canonical quantization is just 
— / £''^x Cc,aa: but this simple relation does not always apply.) Counting a 
factor for each C, or a^, the terms in this interaction go as a^^~'^, a'^'^~^, 
a^'^~'^, a^"'^, and a^"'^, respectively. All these terms are safe for 2u < 4, 
except for the first, which for 2iy > 2 grows exponentially at late times. 

Because of the commutators in Eq. (2), the condition for a safe interac- 
tion is actually less stringent than that it should decay exponentially with 
time, and even a growing term that only involves fields rather than their time 
derivatives, like the first term in Eq. (27), may not destroy the convergence 
of the time integrals. We will now prove the following: 

Theorem The integrals over the time coordinates of interactions converge 
exponentially for t — > oo, essentially as J°° dt/a'^{t) with n > 0, provided 
that in 2v space dimensions, all interactions are of one or the other of two 
types: 

• Safe interactions, that contain a number of factors of a{t) (including 
—2 factors of a for each time derivative and the 2z/ factors of a from 
^/—Deig) strictly less than 21/ — 2, and 
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• Dangerous interactions, which grow at late times no faster than a , 
and contain only fields, not time derivatives of fields. 

These conditions arc evidently met by the interaction (27), irrespective of 
the value of v, and, as we shall see in Section VI, they are satisfied by all 
other interactions in the theories of Section III, but not in all theories. 

Before proceeding to the proof, it should be noted that just as in Eq. (27), 
the space dimensionality 2u enters in the interaction only in a factor ^J—T)etg oc 
a^*^, so the question of whether or not a given theory satisfies the conditions 
of this theorem does not depend on the value of 2z/. Thus this theorem has 
the corollary: 

Corollary The integrals over the time coordinates of interactions converge 
exponentially for t ^ oo, essentially as dt/a"'{t) with n > 0, provided 
that in 3 space dimensions all interactions are of one or the other of two 
types: 

• Safe interactions, that contain a number of factors of a(t) (including 
—2 factors of a for each time derivative and the 3 factors of a from 
yZ—Detg) strictly less than +1, and 

• Dangerous interactions, which grow at late times no faster than a, and 
contain only fields, not time derivatives of fields. 

Here is the proof. As already mentioned, the reason that dangerous 
interactions are not necessarily fatal has to do with how they enter into 
commutators in Eq. (2). Because of the time-ordering in Eq. (2), any failure 
of convergence of the time integrals for t +oo in A^th-order perturba- 
tion theory must come from a region of the multi-time region of integra- 
tion in which, for some r, the time arguments tr, tr+i, ...t]sr, all go to 
infinity together. We will therefore have to count the number of factors of 
a(tr), a(tr+i), • • • a(iAf), treating them all as being of the same order of mag- 
nitude. (This does not take proper account of factors of log a, but as long as 
the integral over tr, tr+i, . . .tN involves a negative total number of factors 
of a, it converges exponentially fast no matter how many factors of log a arise 
from subintegrations.) Now, at least one of the fields or field time derivatives 
in each term in H(ts) with r < s < N must appear in a commutator with one 
of the fields in some other Hjitg') with s < s' < N . So we need to consider 
the commutators of fields at times which may be unequal, but are both late. 
In the sense described above, treating all a{tr), a{tr+i), ...a{tN) as being 
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of the same order of magnitude, if a{t) increases more-or-less exponentially, 
then the commutator of two fields or a field and a field time-derivative goes 
as a~'^'^, while the commutator of two field time-derivatives goes as a~^'^~'^. 

For instance, the unequal-time commutators of the interaction-picture 
fields (14)-(16) are 

C(x, t), C(x', t')] = I cf^p e^P-(---') (Cp WC;(i') - Cp{t'%;{t)) , (28) 

(29) 

'an{^,t),am{^,t')\=5nmj d'" p e'^<^-^'\ap{t)a;it') - ap{t')a;{t)) , 

(30) 

where Ilijki{p) = J2\ ^ij{Pi ^)^kl{Pj -^)- The two asymptotic expansions given 
in Eqs.(21-(23) for each of the fields are both real aside from over-all factors, 
so neither by itself contributes to the commutators. On the other hand, the 
constants CpCp*; ^plp* ^ ^'^d Epa°* are in general complex. (For instance, 
in a strictly exponential expansion, inflation, the phase of CpC,°* is given by 
a factor —e~'^^'^.) The asymptotic expansions of the commutators at late 
times are therefore 



C(xi,ti),C(x2,t2) 



2i I d^^p Im CpC°* 



rt2 



dt' 



t, a^-{t')e{t') 
°° dt[ 



■oo 
oo 



7ij(xi,ii),7fei(x2,i2) 



2i / d^'plii^jkmi^ DpYp 
+p 



„ip-(xi-X2) 



*2 dt' 



2 1 dt I ^2.-2 (^//)^^// / dt 



a2^(t') 



a2-(t;)4 a^^it',)Jt'^ 



0"n(xi,il),Cr^(x2,i2) '^iSnm j d^'^plui 



EpGp* 



„ip-(xi-X2) 



/ 



*2 dt' 
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*2 dt' 

dt'i 



poo fjfll 

^ ' Jt" a2-(t") 

+p2 r r a'-^-\t") dt" 

^ Jt, a2-(ti) Jt, a2-(i'2) A; ^ ' 



+ 



We see that the commutator of two fields vanishes essentially as a"'^'^ for 
late times, and the same is true for the commutator of a field and its time 
derivative, but the commutators of two time derivatives arise only from the 
third terms in the expansions (31)-(33), and therefore go as a"^'^"^. That 



IS, 



C(xi,il),C(x2,t2) 



2i I d'^''pIm[CpC^* 

2 

p 



„ip-(xi-X2) 



f'\'^''-\t')e{t')dt' + 
Ju 



a2'^(ii)e(ii)a2'^(i2)e(i2)4 



and likewise for and (7„. 

Let's now add up the total number of factors of a{tr), a(t,.+i), . . . and a,{ti^] 
in the integrand of Eq. (2), for some selection of terms in the interactions 
H{ts) with r < s < N. Suppose that the selected term in H{ts) contains an 
explicit factor a{ts)^^, and Bg factors of field time derivatives. Suppose also 
that in the inner N — r + 1 commutators in Eq. (2) there appear C com- 
mutators of fields with each other, C commutators of fields with field time 
derivatives, and C" commutators of field time derivatives with each other. 
The number of field time derivatives that are not in these commutators is 

Bs-C- 2C", and these contribute a total -2Y,sBs + 2C" + 4C" factors 
of a. (All sums over s here run from r to N .) In addition, there are 
factors of a that appear explicitly in the interactions, and as we have seen, 
the commutators contribute —2vC — 2vC' — {2v + 2)C" factors of a. Hence 
the total number of factors of a(tr), a(ir-+i)) • • • s^nd a(fjv) in the integrand 
of Eq. (2) is 

# = ^(^-2Ss)-2i/C-(2i/-2)(C' + C") = ^(^-25s-2i/ + 2)-2C , 

s s 

(34) 

in which we have used the fact that the total number C + C" + C" of com- 
mutators of the interactions H{tr), H{tr+i), ■ ■ ■ and H{tN) with each other 
and with the field product Q equals the number of these interactions. Under 
the assumptions of this theorem, all interactions have As — 2Bs < 2v — 2. 



(33) 
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If any of them are safe in the sense that As — 2Bs < 2v — 2, then ^ < 0, 
and the integral over time converges exponentially fast. On the other hand, 
if all of them have Ag — 2Bs = 2v — 2, then under the assumptions of this 
theorem they all involve only fields, not field time derivatives, so the same 
is true of the commutators of these interactions. In this case C > and 
# = — 2C < 0, so again the integral over time converges exponentially fast. 

In counting powers of a, we have held the wave numbers p associated 
with internal lines fixed, like the external wave numbers, because we are 
integrating over time coordinates before we integrate over the internal wave 
numbers. The integrals over time receive little contribution from values of 
the conformal time r = — dt/a satisfying —pr 3> 1 and —qr S> 1, be- 
cause of the rapid oscillation of the integrand, and for theories satisfying the 
conditions of our theorem they also receive little contribution from values of 
r with —pr <C 1 and —qr <C 1, because of the damping provided by negative 
powers of a. (Note that when a{t) increases more or less exponentially, r is 
of the order of —l/aH.) Thus for these theories, we expect the integrals to 
be dominated by times for which — 1/r is in the range from the qs to the pss. 
The question then is whether the integrals over the internal wave numbers 
p are dominated by values of the order of the external wave numbers q? If 
they are, then the results depend only on the history of inflation around the 
time of horizon exit, —qr ~ 1, or in other words, q/a ~ H. 

Any integral over the internal wave numbers will in general take the form 
of a polynomial in the external wave numbers, with coefficients that may be 
divergent, plus a finite term given by a convergent integral dominated by 
internal wave numbers of the same order of magnitude as the fixed external 
wave numbers. An example of this decomposition is given in Section VII. In 
particular, the integral over the wave number associated with an internal line 
that begins and ends at the same vertex does not involve the external wave 
numbers, so its contribution is purely a polynomial in the wave numbers of 
the other lines attached to the same vertex. 

Just as in dealing with ultraviolet divergences in flat space quantum 
field theory, renormalization removes some of these ultraviolet divergent 
polynomial terms, and others are removed by appropriate redefinitions of 
the field operators. (Some examples are given in the next section.) Where 
redefinition of the field operators is necessary, it is only products of the 
redefined "renormalized" field operators whose expectation values may be 
expected to give results that converge at late times. If, after all such renor- 
malizations and redefinitions, there remained ultraviolet divergences in the 
integrals over internal wave numbers, we could conclude that the approxi- 
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mation of extending the time integrals to +00 is not valid, and that these 
integrals can be taken only to some time t late in inflation. The decrease 
of the integrand at wave numbers p much greater than — l/r(t) would then 
provide the ultraviolet cut off that is still needed, but the correlation func- 
tions would exhibit the sort of time dependence that has been found in other 
contexts by Woodard and his collaborators,^ and we would not be able to 
draw conclusions about correlations actually measured at times much closer 
to the present. The possible presence of such ultraviolet divergences that 
are not removed by renormalization and field redefinition is an important 
issue, which merits further study.^^^ But even if such ultraviolet divergences 
are present, it would still be possible to calculate the non-polynomial part 
of the integrals over internal momenta which is not ultraviolet divergent (at 
least in one loop order) even when the time t is taken to infinity. Such a 
calculation will be presented in Section VII. 



V. AN EXAMPLE: EXPONENTIAL EXPANSION 

To clarify the issues discussed at the end of the previous section, we 
will examine a simple unphysical model, along with a revealing class of 
generalizations. 

First, consider a single real scalar field (/^(x, t) in a fixed de Sitter metric. 
In order to implement dimensional regularization, we work in 2z/ space di- 
mensions, letting 1/ ^ 3/2 at the end of our calculation. The Lagrangian 
density is taken as 



2u 2u-2 



(35) 



tt^Many theories are afflicted with infrared divergences, even when t is held fixed. The 
infrared divergences are attributed to the imposition of the unreahstic initial condition, 
that at early times all of infinite space is occupied by a Bunch-Davies vacuum. The 
infrared divergence can be eliminated either by taking space to be finite* or by changing 
the vacuum.® In any case, it is the appearance of uncancelled ultraviolet rather than 
infrared divergences when we integrate over internal wave numbers after taking the limit 
t — » 00 that shows the impropriety of this interchange of limit and integral, because factors 
of 1/a in the integrand are typically accompanied with factors of internal wave numbers, 
so that the 1/a factors do not suppress the integrand for large values of a if the integral 
receives contributions from arbitrarily large values of the internal wave number. 

*This model, and much of the analysis, was suggested to me by R. Woodard, private 
communication. 
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where a oc e with H constant. (This of course can be rewritten as a free 
field theory, but it is instructive nonetheless, and will be generalized later 
in this section to interacting theories.) We follow the usual procedure of 
defining a canonical conjugate field tt = dC/d(p, constructing the Hamilto- 
nian density H = inp — C with ip expressed in terms of tt, dividing H into 
a quadratic part Ho and interaction part Hi, and then replacing tt in Hi 
with the interaction-picture ttj given by tp = [dHo/d-iT]^=^j. This gives an 
interaction 



(An anticommutator is needed in the first term to satisfy the requirement 
that Hi be Hermitian.) This interaction satifies the conditions of the theo- 
rem proved in the previous section for any value of the space dimensionality 
2i': the first term in the square brackets contains — 4 factors of a (count- 
ing a factor a"^ for each time derivative, so it is safe, while the second term 
contains 2u — 2 factors of a, and is therefore dangerous, but it only involves 
fields (including space derivatives), not their time derivatives, so though 
dangerous it still satisfies the conditions of our theorem. 

To first order in A, the expectation value {ip{x,t) ip{x' ,t)) is given by a 
one-loop diagram, in which a scalar field line is emitted and absorbed at the 
same vertex, with the two external lines also attached to this vertex. This 
expectation value receives contributions of three kinds: 

i Terms in which no time derivatives act on the internal lines. This contribu- 

tion is the same as would be obtained by adding effective interactions 
proportional to a^''~^(V(^)^, a^'^'^c/?^, or a?^^'^, all of which satisfy 
the conditions of the theorem of the previous section. Thus it can- 
not affect the conclusion that the integral over the time argument of 
Hi(ti) converges exponentially at ti = -|-oo, so that ((^(x, i) (^(x', t)) 
approaches a finite limit for t — > oo. 

ii Terms in which time derivatives act on both ends of the internal line. This 

produces an effective interaction proportional to c^^ip^ , which violates 
the conditions of our theorem, but it can be removed by adding an 
K^p^ \J —Y)&\.g counterterm in the Lagrangian. (This cancellation is not 
automatic, because the condition of minimal coupling is not enforced 
by any symmetry.) 




(36) 
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iii Terms in which a time derivative acts on just one end of the internal 
line. This produces an effective interaction proportional to o?''^Lpip, 
which violates the conditions of our theorem, and cannot be removed 
by adding a generally covariant counterterm to the Lagrangian. 

To see in detail what trouble is caused by the third type of contribution, 
note that the interaction picture scalar field is given by a Fourier decompo- 
sition like Eq. (16), with coefficient functions'''''' 



gi7r(2i/+l)/4^i/-l/2 



where r is the conformal time 

dt' 



t ait') a{t)H 



(37) 



(38) 



The contribution of the third kind to the expectation value then has the 
Fourier transform 



327r2 



x47r f dtia'^{h)(^ ^_pr,YHi')^_pr, 

Jo P J -co \ati 

{{-qnrHi^^ i-qn )) ' ((-gr)^i7 {-qr 



2k\ 

q J 



Av 



xlm 



*2 



(39) 



Let's see what happens if we evaluate this by integrating first over p and 
then over ti from — oo to late times, or vice versa. 

To integrate first over p, we can change the variable of integration from 

p to z = —pTi, in which case the first derivative with respect to ti can be 
replaced with d/d,ti = {z/aiTi){d/dz) = —Hz{d/dz), while dp/p = dz/z. 
Dimensional rcgularization (with 2ij < 1) makes the function z'^Hi^\z) 
vanish at z ^ oo, while for > it takes the value 2'^r(z/)/7r for z ^ 0, so 



°° , d 
dz — 

dz 



**Here and below we will not be careful to extend factors like A-k to 2v space dimensions. 
This only affects constant terms that accompany any (2z/ — 3)~^ poles. 
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and therefore 



t2u-1 



327r2 



27r\ 
Q. J 



/t ^ 
dti a^''(ti)Im— 
-oo dti 



*2 



(4 



For ti — +00 and t — +00 (that is, r — and ti — 0), the integrand of 
the integral over ti on the second hne has the constant hmit 

4T{ufq 



{{-inrHi'K-in)) {(-qT)-Hm(-qT 



*2 



Thus for t ^ 00, the correlation function (39) does not approach a constant, 
but instead goes as 

/ d-^x e^'>-(--')(^(x,t)^(x',t))i„ ^ ■ (42) 

There is no pole here that prevents continuation to space dimensionality 
2v = 3. Prom this point of view, integrating first over p, the failure of the 
correlation function to approach a finite limit at late times is due to the fact 
already noted, that the integral over p produces an effective interaction that 
does not satisfy the conditions of our theorem. 

But suppose we first integrate over ti from —00 to +00. Now there is no 
problem with convergence at late times, because the original interaction does 
satisfy the conditions of our theorem, but instead we now have a problem 
with the convergence of the integral over p. It will be helpful to divide the 
integral over p into an integral from to Ag, where A ^ 1, and an integral 
from Aq to infinity. The first integral obviously has no ultraviolet divergence, 
and the vanishing of the first time derivative in Eq. (39) for p ^ prevents 
any infrared divergence. In the second integral p and — 1/t are the only 
magnitudes in the problem with which q can be compared, so for t — +00 
and hence r ^ wc can evaluate the correlation function by letting q ^ 
and keeping only the leading term in q. Here again we can use the limiting 
formula (41), now for q ^ instead of r ^ and ti — > 0. The integral over 
ti is then trivial, and we find that for g <C l/r the correlation function is 

/ d^^x e^'i-(— ')(^(x,t)^(x',t))„, ^ ^yj^^o-L'^L £ 7 + • 
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The ultraviolet divergent integral over p is the price we pay for the naugh- 
tiness of taking the limit t — oo before we integrate over p. 

In this model it is clear how to remedy the difficulty of calculating corre- 
lation functions at late times. As already mentioned, the original Lagrangian 
density (35) actually describes a free field theory. This is made manifest by 
defining a new scalar field 



= J ^/l + ^dip , (44) 



for which the Lagrangian density takes the form 

£ = _ 1 v^Zd^ g^^- d^,(pd^(p. (45) 

There is no problem in taking the late-time limit of the correlation function 
j^2r.gjq.(x-x')^(^(x,t)(p(x',t)) — it is just 2'^'' H'^''-'^V{vf /^iir'^q^'' . Prom 
this point of view, the growth of the correlation function (42) at late times is 
a result of our perversity in calculating the correlation function of instead 
of (p. 

Can we find fields whose correlation functions have a constant limit at 
late times in theories that satisfy the conditions of our theorem but arc not 
equivalent to free field theories? The general answer is not known, but here 
is a class of interacting field theories for which such "renormalized" fields 
can be found. This time we consider an arbitrary number of real scalar fields 
(pn (x, t) in a fixed de Sitter metric. The Lagrangian density is taken to have 
the form of a non-linear cr-model: 

^ = -\Y^ V-Detg g^"' {Snm + ><Knm{v)) d/^ipn dt,(pm , (46) 



nm 



where Knm{'p) is an arbitrary real symmetric matrix function of the ipn\ A 
is a coupling constant; and again a cx e^* with H constant. The Hamilto- 
nian derived from this Lagrangian density docs satisfy the conditions of the 
theorem of Section IV, whatever the function Knm{(p). 

To first order in A, the same problem discussed earlier in this section 
arises from graphs in which an internal line of the field (fn is emitted and 
absorbed from the same vertex, with a time derivative acting on just one 
end of this line. Depending on what correlation function is being calculated, 
the contribution of such graphs is proportional to various contractions of 
partial derivatives of the function 
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Suppose we make a redefinition of tlie fields of first order in A: 



<f>n = Vn- AA„((^) . (48) 

This changes the matrix K to 

J^nmVP) = J^-nmif) H ^ 1 ^ , (49) 

and so 

Thus the fields are renormalized, in the sense that to first order in A 
correlation functions have finite limits at late times, provided that 

This can be solved by first solving the Poisson equation 

^ a^^(y) _ 1 ^ dAM 

and then solving a second Poisson equation 

2:fM»:)=_^„(,)_fM. (53) 

^ d(pnd(pn difm 

Thus for at least to first order in this class of theories, it is always possible 
to find a suitable set of renormalized fields. 

Because we can take the limit t ^ oo only for the correlation functions 
of suitably defined fields (such as (pn in our example), the question naturally 
arises, whether these are the fields whose correlation functions we want to 
calculate. The answer is conditioned by the fact that astronomical observa- 
tions of the cosmic microwave background or large scale structure are made 
following a long era that has intervened since the end of inflation, during 
which things happened about which we know almost nothing, such as re- 
heating, baryon and lepton synthesis, and dark matter decoupling. The only 
thing that allows us to use observations to learn about inflation is that some 
quantities were conserved during this era, while fluctuation wave lengths 
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were outside the horizon. These are the only quantities whose correlation 
functions at the end of inflation can be interpreted in terms of current ob- 
servations. In the classical limit, the quantities that are conserved outside 
the horizon arc and 7^ , but we don't know whether this will be true when 
quantum effects are taken into account. Still, we can expect that quantities 
are conserved only when there is some symmetry principle that makes them 
conserved, and whatever symmetry principle keeps some quantity conserved 
from the end of inflation to the time of horizon re-entry is likely also to keep 
it conserved from the time of horizon exit to the end of inflation. So we may 
guess that the quantities whose correlation functions we will need to know 
are just those whose correlation functions approach constant limits at the 
end of inflation. 



VI. DANGEROUS INTERACTIONS IN INFLATIONARY 

THEORIES 

We now return to the semi-realistic theories described in Section III. We 
will show in this section that all interactions are of the type called for in the 
theorem of Section IV; that is, they are all safe interactions that (in three 
space dimensions) do not grow exponentially at late times (and in fact are 
suppressed at late times at least by a factor a^^), or dangerous interactions 
containing only fields and not their time derivatives, which grow no faster 
that a at late times. Fortunately, as noticed by Maldacena^ in a different 
context, for this purpose it is not necessary to solve the constraint equations 
(7) and (8), which are quite complicated especially when the an fields are 
included. Inspection of these equations shows that when we count ^5 iij, 
and an as of order a~^, the auxiliary fields N — 1 and iV* arc both also of 
order a^^.-tW This is apparent in the first-order solution (9) of the constraint 
equations, but it holds to all orders in the fields. To calculate the quantity 
E^iE'j - {E\f in Eq. (5), we note that 



dt 



_ Vv^iV^. + V.iV*) . (54) 



The first term HS^'j is of order zero in a, while all other terms are of order 
a~^, so 

E^iE'j - {E\f = -6ij2 - \2HC, - AHVkN'' + 0(0"^) (55) 

In counting powers of a, note that the three-dimensional affine connection and Ricci 
tensor are independent of a, so the curvature scalar 7?*^' goes as a~^. For instance, for 
■yij = 0, we have R^^^ = -a-^e-^^iW^C + 2(VC)^). 
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= in = 0.) The terms 
as a consequence of the constraint 



(In deriving this result, we note that 
in (5) of first order in iV — 1 ah cance. 
equation (8), while terms of second order in AT — 1 in Eq. (5) (and in partic- 
ular in a^e^'^ip^ /2N and —ZH^a^e^^/2N) are suppressed by at least a factor 
a^(a^^)^, and arc therefore safe. Therefore we can isolate all terms that are 
potentially dangerous by setting = 1, and find 



+ (^^-a ^^[exp(-7)]*-'^5iC7„0j(T„ 



+ 0(a-^), 



(56) 



We note that e^^V^A^*' = dkie^^N''), so this term vanishes when integrated 
over three-space, and therefore makes no contribution to the action. The 
term proportional to <^ can be written 

-6a^e^^HC = ^ (-2a^ e^^) + a^e^^(6H^ + 2H) . 

The first term vanishes when integrated over time, so it gives no contribution 
to the action. To evaluate the remaining terms we use the unperturbed 
infiaton field equation, which (with SttG = 1) gives H = —(p /2, and the 
Priedmann equation, which gives 6iJ^ = 2V+i^'^ . We then find a cancellation 

-V - 3if2 + + 6if^ + 2H = Q . 

Aside from terms that make no contribution to the action, the Lagrangian 
density is then 



=3C 



'-2e-2C[exp(-7)]*^X^^^n5,^„ 



+ 0(a 



(57) 



We see that, at least in this class of theories, the dangerous terms that are 
not suppressed by a factor grow at most like a at late times, and involve 
only fields, not their time derivatives, as assumed in the theorem of Section 
III. 

It remains to be seen if in these theories, after integrating over times and 
taking the limit t — ^ oo, the remaining integrals over internal wave numbers 
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are made convergent by the same counterterms that ehminate ultraviolet di- 
vergences in fiat spacetime, and if not, whether they can be made convergent 
by suitable redefinitions of the fields ( and 'jij appearing in the correlation 
functions. This is left as a problem for further work. 

Not all theories satisfy the conditions of the theorem of Section IV. For 
instance, a non-derivative interaction ^/—DetgF{a) of the a fields would 
have +3 factors of a, and hence would violate the condition that the total 
number of factors of a (counting each time derivative as -2 factors) must 
be no greater than +1. The a fields must be the Goldstone bosons of some 
broken global symmetry in order to satisfy the conditions of our theorem in 
a natural way. 

VII. A SAMPLE CALCULATION 

As an application of the formalism described in this paper, we will now 
calculate the one-loop contribution to the correlation function of two ^ fields, 
which is measured in the spectrum of anisotropies of the cosmic microwave 
background. As already mentioned, in the class of theories described in 
Section III, this two-point function is dominated by a matter loop, because 
there are many types of matter field and only one gravitational field. We 
first consider the contribution of second order in the interaction (27). It 
saves a great deal of work if we use the interaction-picture field equations 
(11) and (13) to put this interaction in the form 

Hcaait) = -Jd^x >Ccaa(x, t) = A{t) + B{t) (58) 

where 

A = -2eHa^ ^ / d^x al V^^C (59) 

n •' 

B = lljd'x (^-ea^V-^c) (^(Va„)^ + ^aV^) . (60) 

In general, for an interaction Hamiltonian of the form (58), Eq. (2) can be 
put in the form 

(Qit)) = dtN-1 ■■■ dh 

»r n J—OO J — OO J—OO 



N=0 
X 



{[Hi{ti),[Hi{t2),---[HiitN),Q\t)] ■■■]]) , (61) 
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where 



- JB{t) 



-iB{t) _ 



A{t)+i[Bit),A{t)]+-m),B{t)]+... 
(62) 



Q\t)e 



-iBit) ^ (^I 



Q\t)+i[B{t),Q'{t)]--[B{t),[B{t),Q'{t)]]+--- ■ 

(63) 



To second order in an interaction of the form (58), the expectation value is 
then 



m)h 



dt2 dh(\Aiti),\Ait2),Q\t) 

-OO J — OD ^ '- '- 

/* dh (\\B{h),A{h) + B{h)/2\,Q'{t) 
B{t),\Bit),Q'{t) 



(64) 



The Fourier transform of the second-order term in the expectation value of 
a product of two <^s is then 



vac, m 



J d^x e^'i-(''-''')(vac, in|c(x, t) C(x', t) 

= -^^^Re/* a''it2)e{t2)Hit2)dt2 
Q J—oo 

X a^{ti)e{ti)H{ti)dti 

J—oo 

x4(ii)C(t) (4(i2)c*(t) -ut)c;{t2)) 

xM jd^pj d^p' (^3(p + p' + q) 

X (t 1 ) (T* {t2 ) Opi {ti ) &*, {t2 ) 



(2vr) 



-M J d^p J d^p' 53 (p + p' + q) 



4g4 

x(p.p')2|ap(i)nap,(i)l' 
+ ... 



(65) 



where is the number of a fields. We have shown here explicitly the 
contribution of the first and third lines on the right-hand side of Eq. (64). 
The dots represent one-loop contributions of the second line, in which [B, A+ 
B/2\ plays the role of a CC"''^ "seagull" interaction, as well as one-loop terms 
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of first order in the CC*^"^ terms in Eq. (5), in both of which the integral 
over internal wave number is q-independent, plus counterterms arising in 
first order from interactions that cancel ultraviolet divergences in flat space, 
including ^/—Detg R^y and ^J—T)e\,g B? terms in the Lagrangian density 
that arc not included in Eq. (5). 

Though it has not been made explicit in this section, we use dimensional 
regularization to remove infinities in the integrals over p and p' at interme- 
diate stages in the calculation, and we now assume that the singularity as 
the number of space dimensions approaches three is cancelled by the terms 
in Eq. (65) represented by dots, leaving it to future work to show that this 
is the case. Then these integrals are dominated hy p p' k. q. As we have 
seen, the integrals over time are then dominated by the time tq of horizon 
exit, when q/a{tq) ~ H{tq). For simplicity, we will assume (for the first time 
in this paper) that the unperturbed infiaton field if{t) is rolling very slowly 
down the potential at time tq, so that the expansion near this time can be 
approximated as strictly exponential, a{t) oc e^*. Then the wave functions 



are 



aq{t)c,a°e-"i^[l + iqT) , 
Cq{t)cCqe-''^{l + iqT), 



where r is the conformal time 



r 



dt 



It ait) ' 

and the wave functions outside the horizon have modulus 



H'jtq) 

2(27r)3 g3 ' 



IC 



q\ 



H\tq 



2(27r)3e(i,)g3 



Using these wave functions in Eq. (65) gives 



'3xe^'i(^-^')(vac,m 



vac, m 



(27r)2 



pp 



C(x,t)C(x',t) 

j d^p j dV<5^(p + p' + q) 
(p-pO' 



+ 



q'^{p+p' + q) 16g^p3p'3 



+ 



(66) 



with the dots having the same meaning as in Eq. (65). 
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Simple dimensional analysis tells us that when the integral over internal 
wave numbers of the first term in square brackets is made finite by dimen- 
sional regularization, it is converted to 



/ j dy + p' + q) ^ g^+^F(<5) , 

J J p + p +q 



(67) 



where 5 is a measure of the difference between the space dimensionality 
and three. The ultraviolet divergences in this integrals for 5 = gives the 
function F{5) a singularities as 5 — 0: 



Fa 



(68) 



so that in the limit 6 = 



j d^p j d^p' (53(p + p' + q) 



pp 



p+p' + q 



Folnq + L 



(69) 



where L is a divergent constant. We can easily calculate the coefficient Fq 
of the logarithm. For this purpose, we note that, in general, 

/" 3 /" 3 / 3 / / 2tT f°° fP+1 I I , 

d p d p 6 {p + p +q}f{p,p ,q) = — / pdp pdpf{p,p,q) 
J J q Jo J\p-i\ 

(70) 

To eliminate the divergence in the integral over p and p', we multiply by q 
and differentiate six times with respect to q. A tedious but straightforward 
calculation gives 



d^ 
dq^ 



Qjd^pj d^p' 6^{p + p' + <i] 



PP 



p + p' + q 



Stt 
Q 



Comparing this with the result of applying the same operation to Eq. (69) 
then gives Fq = — 7r/15. 

In contrast, the integral of the second term in square brackets in Eq. (66) 
is a sum of powers of q with divergent coefficients, but with no logarithmic 
singularity in q. (This term would be eliminated if we calculated the ex- 
pectation value of a product of fields C = exp(— zi?)C cxp(ii3) instead of (.) 
The terms represented by dots in Eq. (65) make contributions that are also 
just a sum of powers of q with divergent coefficients. We are assuming that 
all ultraviolet divergences cancel, but we cannot find resulting finite power 
terms without knowing the renormalized coefficients of the ^/ —'DetgW^ R^xv 
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and ^/—DetgR^ terms in the Lagrangian density. So we are left with the 
result (now restoring a suitable power of SttG) that 



/ 



Tr(8TrGH^{tq)yM 



vac, in 

2 



15(27r)3g3 



C(x,t)C(x',t) 



Ing + C 



vac, in 



(71) 



with C an unknown constant. This may be compared with the classical (and 
classic) result, that in slow roll inflation this correlation function takes the 
form 



^3^giq.(x-x') 



vac, m 



C(x,t)C(x',t) 



vac, m 



87rGH^{tq) 
4(27r)3|e(i,)|Q3 



(72) 



The one-loop correction (71) is smaller by a factor of order 8'KGH'^M\e(tq)\, 



so even if Af is 10^ or 10^ this correction is likely to remain unobservable. 
Still, it is interesting that even in the extreme slow roll limit, where H{tq) 
and e{tq) are nearly constant, the factor In q gives it a different dependence 
on the wave number q. 
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APPENDIX: THE IN-IN FORMALISM 



1. Time Dependence 

First, it is necessary to be precise about the origin of the time-dependence 
of the fluctuation Hamiltonian in applications such as those encountered in 
cosmology. Consider a general Hamiltonian system, with canonical variables 
(/>a(x, i) and conjugates 7ra(x, t) satisfying the commutation relations 



(j)a{yi,t),TTb{y, t) 



and the equations of motion 

4>a{^,t)=i\H[(t>{t),TT{t)\,(t>a{^,t) 



baix, t),(f)b{y,t) 



7ra(x, t) 



7ra(x,t),7r6(y,t) 
(A.l) 



0, 



H[(P{t),7r{t)],Tra{x,t) . 

(A.2) 

Here a is a compound index labeling particular fields and their spin com- 
ponents. The Hamiltonian H is a functional of the (pai^^t) and 7ra(x, t) at 
fixed time t, which according to Eq. (A.2) is of course independent of the 
time at which these variables are evaluated. 

We assume the existence of a time-dependent c-number solution ^a(x, t), 
7fa(x, t), satisfying the classical equations of motion: 



7ra(x, t) 



5(j)a{yi,t) 



(A.3) 



(57ra(x,t) 

and we expand around this solution, writing 

(?!)„(X, t) = 4(X, t) + S(f)a{x, t) , 7r„(x, t) = 7fa(x, t) + STTaix, t) . (A.4) 

(In cosmology, 0a would describe the Robertson-Walker metric and the 
expectation values of various scalar fields.) Of course, since c- numbers com- 
mute with everything, the fluctuations satisfy the same commutation rules 
(A.l) as the total variables: 



^a(x, t),5nb{y,t) 



^a(x, i),(506(x, t) 



67ra{-^,t),6TTb{x,t) 

(A.5) 

When the Hamiltonian is expanded in powers of the perturbations (5(^a(x, t) 
and 5TTa{x,t) at some definite time t, wc encounter terms of zcroth and first 
order in the perturbations, as well as time-dependent terms of second and 
higher order: 



H[<j>{t),n{t)] = Hm,7rit)] + Y^ 



5H[^it),Ttit)] 



6(l)a{y:,t) 
+ H[S4>{t),S7Tity,t] , 
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where H[6(l){t),d7r{t);t] is the sum of all terms in H[(j){t) + S(l){t),7r{t) + d7r{t)] 
of second and higher order in the 64>{x,t) and/or 6Tr{x,t). 

Now, although H generates the time-dependence of (pai^, t) and 7ra(x, t), 
it is H rather than H that generates the time dependence of S(j)a{x,t) and 
(57ra(x,i). That is, Eq. (A.2) gives 



while Eqs. (A. 5) and (A. 3) give 



TTa{y:,t)+STra(yi,t) = 1 H[(f>{t) ,Tr{t)], STra(p<.,t) 



5H[4>{t),Tt{t)] 



^ . S(f)b{y,t) 
Subtracting, we find 



Sn{y, t) 



(57rb(y, t),S(l)aiyi, t) 



^7r6(y,t),57ra(x, t) 



S<j)a{x,t) = i[H[(Pit),7rity,t],Scl)a{x,t)\ , S7ra{x,t) = i[H[<Pit),n{t)-t],67Taix,t) 

(A.7) 

This then is our prescription for constructing the time-dependent Hamilto- 
nian H that governs the time-dependence of the fluctuations: expand the 
original Hamiltonian H in powers of fluctuations 5(p and 5tt, and throw 
away the terms of zeroth and first order in these fiuctuations. It is this 
construction that gives H an explicit dependence on time. 

2. Operator Formalism for Expectation Values 

We consider a general Hamiltonian system, of the sort described in the 
previous subsection. It follows from Eq. (A.7) that the fiuctuations at time 
t can be expressed in terms of the same operators at some very early time 
to through a unitary transformation 

5(t>a{t) = U-\t,to)5Mto)Uit,to) , STTait) = U-\t,to)6Tra{to) U {t,to) , 

(A.8) 

where U{t,to) is defined by the differential equation 



di 



U{t, to) = -i H[6(l)it), 6Tr{ty,t] U{t, to) 



and the initial condition 



Uito,to) 



(A.9) 
(A.10) 
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In the application that concerns us in cosmology, we can take to = —oo, hy 
which we mean any time early enough so that the wavelengths of interest 
are deep inside the horizon. 

To calculate U{t, to), wc now further decompose H into a kinematic term 
Hq that is quadratic in the fluctuations, and an interaction term Hj: 

H[S(j){t),6Tr{ty,t] = Ho[S(j)it),6Tr{ty,t]+Hi[6(j){t),6Trit);t] , (A.ll) 

and we seek to calculate U as a power series in Hj. To this end, we intro- 
duce an "interaction picture": we define fluctuation operators S(/)j^(t) and 
S7r^{t) whose time dependence is generated by the quadratic part of the 
Hamiltonian: 

5<pi{t)=i[Ho[Scp\t),67r'{ty,t],6<piit)] , 67ri{t) = i[Ho[5^\t),67r' {ty,t],67riit) 

(A.12) 

and the initial conditions 

S^iito) = SMto) , Siriito) = SiTaito) . (A.13) 

Because Hq is quadratic, the interaction picture operators are free fields, 
satisfying linear wave equations. 

It follows from Eq. (A.12) that in evaluating Ho[d(f)^ ,67:^ -,1] we can take 
the time argument of Scj)^ and Stt^ to have any value, and in particular we 
can take it as to, so that 

Ho[S(f>\t),S7r^{tyt] = Ho[S(l){to),S7T{toyt] , (A.14) 

but the intrinsic time-dependence of Hq still remains. The solution of 
Eq. (A.12) can again be written as a unitary transformation: 

S4>i{t) = UoHt, to)SMto)Uo{t, to) , Sniit) = UQ\t, to)S7raito)Uoit, to) , 

(A.15) 

with Uo defined by the differential equation 

jUo{tM) = -iHo[S(j>{toySTr{toy,t] C/o(t,<o) (A.16) 
and the initial condition 

C/o(to,to) = l- (A.17) 
Then from Eqs. (A. 9) and (A.16) we have 
d 



dt 



Uo\t,to)U{t,to) =-iUo\t,to)Hi[d(f>{to),S7r{toyt]U{t,to) . 
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Using Eq. (A. 15), this gives 



U{t,to) = Uo{t,to)F{t,to) , (A.18) 

where 

-F{t, to) = -iHi{t)F{t, to) , F{to, to) = l. (A.19) 

and Hj{t) is the interaction Hamiltonian in the interaction picture: 

Hi{t) = Uo{t,to)Hi[5(f){to),6Tr{to);t]Uo\t,to) = Hi[6cl)\t),67T^ {t);t] 

(A.20) 

The solution of equations like (A.19) is well known 

F{t,to) = Te^p(^-i Hj{t)dt^ (A.21) 

where T indicates that the products of Hjs in the power series expansion of 
the exponential are to be time-ordered; that is, they are to be written from 
left to right in the decreasing order of time arguments. The solution for the 
fluctuations in terms of the free fields of the interaction picture is then given 
by Eqs. (A.8) and (A.15) as 



Q{t) = F-^{t,to)Q'it)Fit,to) 



Texi)(^i J^^ Hi{t)dt^ Q\t) T ex^ {-i Hi{t) dt^ 



A-22) 



where Q{t) is any (5(/)(x, t) or (y7r(x, t) or any product of the 5(j)s and/or Jtts, 
all at the same time t but in general with different space coordinates, and 
Q^{t) is the same product of 6(j)\^,t) and/or (57r^(x, t). Also, T denotes 
anti-time-ordering: products of Hjs in the power series expansion of the 
exponential are to be written from left to right in the increasing order of 
time arguments. 

3. Diagrammatic Formalism for Expectation Values 

We want to use Eq. (A. 22) to calculate the expectation value {Q{t))ol 
the product Q{t) in a "Bunch-Davies" vacuum, annihilated by the positive- 
frequency part of the interaction picture fluctuations d(p^ and Stt^ . We can 
use the familiar Wick theorem to express the vacuum expectation value of 
the right-hand side of Eq. (A. 22) as a sum over pairings of the 6(p^ and 6tt^ 
with each other. (This of course is the same as supposing the interaction- 
picture fields in Hi{t) and Q^{t) to be governed by a Gaussian probability 
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distribution, except that the order of operators in bihnear averages has to 
be the same as the order in which they appear in Eq. (A. 22).) Expand- 
ing Eq. (A. 22) as a sum of products of bihnear products leads to a set of 
diagrammatic rules, but one that is rather complicated. 

In calculating the term in [Q) of A^^th order in the interaction, wc draw 
all diagrams with N vertices. Just as for ordinary Feynman diagrams, each 
vertex is labeled with a space and time coordinate, and has lines attached 
corresponding to the fields in the interaction. There are also external lines, 
one for each field operator in the product Q, labeled with the different 
space coordinates and the common time t in the arguments of these fields. 
All external lines are connected to vertices or other external lines, and all 
remaining lines attached to vertices are attached to other vertices. But there 
are significant differences between the rules following from Eq. (A. 22) and 
the usual Feynman rules: 

• We have to distinguish between "right" and "left" vertices, arising 
respectively from the time-ordered product and the anti-time-ordered 
product. A diagram with N vertices contributes a sum over all 2''^ 
ways of choosing each vertex to be a left vertex or a right vertex. 
Each right or left vertex contributes a factor —i or +i, respectively, as 
well as whatever coupling parameters appear in the interaction. 

• A line connecting two right vertices or a right vertex and an external 
line, in which it is associated with field operators A(x, t') and S(y, t"), 
contributes a conventional Feynman propagator {T{A{-x.,t')B{y,t"}). 
(It will be understood here and below, that in calculating propagators 
all fields A, B, etc. are taken in the interaction picture, and can 
be Stp^s and/or Stt^s.) As a special case, if B is associated with an 
external line then t" = t, and since t' < t, this is {B{y,t)A{-x.,t')). 

• A line connecting two left vertices, associated with field operators 
A{x,t') and B{y,t"), contributes a propagator {f{A{yL,t')B{y,f"}). 
As a special case, if B is associated with an external line then t" = t, 
and this is {A{x,t')B{y,t)). 

• A line connecting a left vertex, in which it is associated with a field 
operator A(-x,t'), to a right vertex, in which it is associated with a 
field operator B{y,t"), contributes a propagator {A{yi,t')B{y,t")). 

• We must integrate over all over the times t', t" , . . ., associated with the 
vertices from to to t, as well as over all space coordinates associated 



31 



with the vertices. 



We must say a word about the disconnected parts of diagrams. A vac- 
uum fluctuation subdiagram is one in which each vertex is connected only to 
other vertices, not to external lines. Just as in ordinary quantum field the- 
ories, the sum of all vacuum flTictiiation diagrams contributes a numerical 
factor multiplying the contribution of diagrams in which vacuum fluctua- 
tions are excluded. But unlike the case of ordinary quantum field theory, 
this numerical factor is not a phase factor, but is simply 



Texp i 



Jto 



dt 



Texp 



; fHj{t) 
Jto 



dt 



1 . 



(A.23) 



Hence in the "in-in" formalism all vacuum fluctuation diagrams automati- 
cally cancel. Even so, a diagram may contain disconnected parts which do 
not cancel, such as external lines passing through the diagram without in- 
teracting. Ignoring all disconnected parts gives what in the theory of noise 
is known as the cumulants of expectation values,^*^ from which the full ex- 
pectation values can easily be calculated as a sum of products of cumulants. 

4. Path Integral Derivation of the Diagrammatic Rules. 

It is often preferable use path integration instead of the operator for- 
malism, in order to derive the Feynman rules directly from the Lagrangian 
rather than from the Hamltonian, or to make available a larger range of 
gauge choices, or to go beyond perturbation theory. Going back to Eq. (1), 
and following the same reasoning^^ that leads from the operator formalism 
to the path-integral formalism in the calculation of S-matrix elements, we 
see that the vacuum expectation value of any product Q{t) of 54)S and Jtts 
at the same time t (now taking to = — oo) is 



dSlTRai'^, t') 



I n ds^La{^,t') n ""^'^^T"^ n dS4>naM n 

x,t',a IT ^' , 



X exp 



x,t',a 



K.,t' ,a 



dt' 



J2 j d^x5^La{^,t')5lTLa{^,t') - H[5^L{t'), 5lTL{t')\ t'] 



expjiy* dt! d^^S^Ra{^,t')STrRa{^,t')-H(6M^,t'),^^R{^,t'y,t') | 

'[lS(S(j)Lai^,t) - 5^i?a(x,t)) 6(5'KLai^,t) - 5'KRa{^,t)^ Q 54>L{t),5'KL{t) 

*^[^</.l(-oo)] *o[<50ii(-oo)] . (A.24) 
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Here the functional ^'o['^^] is the wave function of the vacuum, 

*o[0(-oo)] oc exp I -^X! / / d^y^abi^,y)S4'a{^,-oo)6(j)biy,-oi 

= exp ^-i £^dt' e^'' J d^x J d'y£abi^,y)dMt')SMt')^ (A-25) 

where £ab is ^ positive-definite kernel. For instance, for a real scalar field of 
mass m, 

f (x, y) = I d^p e^P-("-y) + m2 . (A.26) 

As is well known, if the Hamiltonian is quadratic in the canonical con- 
jugates SiTa with a field-independent coefficient in the term of second order, 
then we can integrate over the Sna by simply setting 5(pa = dH/dSiTa, and 
the quantity Ea S4>ait')STra{t') - H(S(l){t'), 6Tr{t'); t'^ in Eq. (A.24) then be- 
comes the original Lagrangian. We will not pursue this here, but will rather 
take up a puzzle that at first sight seems to throw doubt on the equivalence 
of the path integral formula (A.24), when we do not integrate out the tts, 
with the operator formalism. 

The puzzle is that, although the propagators for lines connecting left 
vertices to each other or right vertices to each other or left or right vertices 
to external lines are Greens functions of the sort that familiarly emerge 
from path integrals, what are we to make of the propagators arising from 
Eq. (A. 22) for lines connecting left vertices with right vertices? These are not 
Greens functions; that is, they are solutions of homogeneous wave equations, 
not of inhomogeneous wave equations with a delta function source. As 
we shall see, the source of these propagators lies in the delta functions in 
Eq. (A.24). It is these delta functions that tie together the integrals over 
the L variables and over the R variables, so that the expression (A. 18) does 
not factor into a product of these integrals. 

In analyzing the consequences of Eq. (A.24), it is convenient to condense 
our notation yet further, and let a variable ^n(i) stand for all the 6(j)a{x,t) 
and (^7ra(x, i), so that n runs over positions in space and whatever discrete 
indices are used to distinguish different fields, plus a two- valued index that 
distinguishes Scf) from Sir. With this understanding, Eq. (A.24) reads 



t'n " t',n " 
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X exp < —I 



dt' L[^L{t')AL{t'); t')Y^v[i f_^dt' L[^R{t'),iR{t')- t' 

'^^{iLn{t)-iRn{t)) \ Q(CL(i))^o(a(-Oo))*o(Cil(-Oo)) , (A.27) 



where 



(A.28) 

To expand in powers of the interaction, we spht L into a term Lq that is 
quadratic in the fluctuations, plus an interaction term —Hj: 



L = Lq — Hi 



(A.29) 



where 



Lorn'), i{t'); = E / d^xdM^,t')S7raM - Ho{Sct>it'), dirit'); t') . 

a •' 

(A.30) 

As in calculations of the S-matrix, we will include the argument of the 
exponential in the vacuum wave functions along with the quadratic part of 

the Lagrangian, writing 



dt'\Lo[iR{t'),iR{t');t'] 

-oo I 

+ f E / J dSSab{^,y)S(l>R^i-x,t')S(l>Rbiy,t')\ 

ab ) 

^lj2J2^nt',mt"^Rn{t')^Rn'{t") , (A.31) 

nn' t',t" 

* dt'hoiut'), Ut'y, t'] 



f E /^'^ I d^ySabi^,y)6<j>LaM6cj)n>iy,t') 

ab 



= 9 E E '^nt',n't" ^Ln{t') ^Ln'{t") 



(A.32) 



nn' t',t" 



The vacuum wave function is the same for and ^r, but it is combined here 
with an exponential exp(— z J Lq) for the ^Ln an exponential exp(+z J Lq) 
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for the ^Rn, which accounts for the different signs of the ie terms in Eqs. (A. 31) 
and (A. 32). (The factor e^* in Eq. (A. 25) is effectively equal to one for any 
finite t', and has tlicrcforc been dropped.) We also express the product of 
delta functions in Eq. (A. 27) as a Gaussian: 

UK^Lnit) - Unit)) OC exp ( -i i^Lnit) - ?iln(t)) ) 

J212^nt',n't"{Un{t') " Un{t')) {Un'{t") " (i")) jA.33) 
nn' ft" J 



= exp 

where 



Cut', n't" = ^Snn' S{t' - t) 5{t" - t) , (A.34) 

and e' is another positive infinitesimal. 

Following the usual rules for integrating a Gaussian times a polynomial, 
the integral is given by a sum over diagrams as described above, but with 
a line that connects right vertices with each other (or with external lines) 
contributing a factor —i^^fn't"^ * ^^^^ ^^^^ connects left vertices with each 
other (or with external lines) contributing a factor iA^^^,^„, and a line 
that connects a right vertex where it is associated with ^n(t') with a left 
vertex associated with ^n'{t") contributing a factor iA^^^n„, with the As 
determined by the condition 



W'^-C C W -iA^^ iA«^ ^ _ 1 
C -iV^ -C ] \ i{A^^y iA^^ J ~ I 1 



(A.35) 



This must hold whatever tiny value we give to e', and so 

jyH^HR ^ ^ ^ v^A^^ = 1 , (A.36) 

T^Rl^RL ^ Q ^ (a^^) ^ = , (A.37) 

CA^^ = CA^^ , CA^^ = -C(A^^)T . (A.38) 

The first Eq. (A.36) is the usual inhomogeneous wave equation for the 
propagator, whose solution as well known is 

-iA^^^n't" = {T{Ut') in'{t")}) , (A.39) 
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with the time-ordering dictated by the +ie in Eq. (A. 31). The second 
Eq. (A. 36) is the complex conjugate of the first wave equation, whose solu- 
tion is the complex conjugate of Eq. (A.39): 

iA^i;^^,,„ = {f{Ut')U'(.t")}) . (A.40) 

Eqs. (A.39) and (A.40) thus give the same propagators for lines connecting 
right vertices with each other or with external lines, and for lines connecting 
left vertices with each other or with external lines, as we we encountered 
in the operator formalism. Equations (A.37) tell us that A^^ and (A^^)""" 
satisfy the homogeneous versions of the wave equations satisfied by A^'^ and 
A^^, but to find A^^ we also need an initial condition. This is provided by 
the first of Eqs. (A. 38), which in more detail reads 

iA^^it, h) = lA^^^it, = {f{Ut)im{t2)]) = {U{t2)Ut)) , (A.41) 

in which wc have used the fact that t > t2- This, together with the first of 
Eqs. (A.37), tells us that 

iA^^ih,t2) = {U{t2)Uti)) , (A.42) 

which is the same propagator for internal lines connecting right vertices with 
left vertices that we found in the operator formalism. 

5. Tree Graphs and Classical Solutions. 

We will now verify the remark made in Section I, that the usual approach 

to the calculation of non-Gaussian correlations, of solving the classical field 
equations beyond the linear approximation, simply corresponds to the cal- 
culation of tree diagrams in the "in-in" formalism. This is a well-known 
result^^ in the usual applications of quantum field theory, but some modifi- 
cations in the usual argument are needed in the "in-in" formalism, in which 
the vacuum persistence fTinctional is always unity whether or not we add a 
current term to the Lagrangian. 

We begin by introducing a generating functional W[j, t, g] for correlation 
functions of fields at a fixed time t: 



^W[J,t,g]/g ^ / ^^^^ 



e 9 




(A.43) 



where Ja is an arbitrary current, and g a real parameter, with the sub- 
script g indicating that the expectation value is to be calculated using a 
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Lagrangian density multiplied with a factor 1/g. (This is different from the 
usual definition of the effective action, because here we are not introducing 
the current into the Lagrangian.) The quantity of physical interest is of 
course W[J, t, 1], from which expectation values of all products of fields can 
be found by expanding in powers of the current. 

Using Eq. (A. 27), we can calculate W as the path integral 

= y jYldTTL lll^cl>R JUSnR 

X exp ^— i J dt' -L[S(l)L,S7rL;t']^ 
X exp J dt' -L[S(l)R,6TrR;t']J 

X ^v<«:['^<^l(-00)] *vac[(^<^i?(-00)] (A.44) 

The usual power-counting arguments^^ show that the L loop contribution 
to W[J,t,g] has a ^-dependence given by a factor g~^. For g ^ 0, W is 
thus given by the sum of all tree graphs. The integrals over Scpi, 5ttl, Scpi^, 
Sttl are dominated in the limit ^ — >^ by fields where L is stationary, i.e., 
where 

with (J^ciassicai ^^ciassicai ^]^g solutions of the classical field equations 
with the initial conditions that the fields go to free fields such as (14)-(16) 
satisfying the initial conditions (20) at t ^ — oo. Since the L and R fields 
take the same values at this stationary point, the action integrals cancel, 
and we conclude that 

W[J, t, 1]] = ^ / d^x -'-'(x, t) Ja{^) . (A.45) 

J zero loops '-^ J 

Expanding in powers of the current, this shows that in the tree approxima- 
tion the expectation value of any product of fields is to be calculated by 
taking the product of the fields obtained by solving the non-linear classical 
field equations with suitable free-field initial conditions, as was to be proved. 
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